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Abstract. We present the results of fast temporal op- 
tical spectroscopy and photometry of the AM Her type 
system BY Cam obtained in February 1990 and March 
1991. Emission line profiles show a complex structure and 
are strongly variable. The radial velocity studies of a sharp 
component detected in Balmer lines in March confirms the 
non-synchronism of the system. Different possible emit- 
ting regions are discussed to explain the characteristics of 
the multiple line components. It is shown from a computa- 
tion of the line velocities and widths that the 'horizontal' 
stream cannot explain the lines of intermediate widths 
observed in this system. A critical review of the different 
estimations of the periods found in BY Cam is presented. 
Additional periods are revealed from a re-analysis of pre- 
vious optical polarimetric and UV spectroscopic data. We 
show that because of the slight asynchronism between the 
white dwarf rotation and the orbital period, significant 
changes in the accretion geometry introduce a bias in the 
period determination depending on the length of the ob- 
servations. Large phase variations are shown to exist which 
are well reproduced by a phase-drift model in which a 
magnetic dipole rotates in the orbital frame with a period 
of 14.5±1.5 days. 
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1. Introduction 

Polars (AM Her type systems) are cataclysmic binaries 
in which a magnetic white dwarf accretes matter from a 
red dwarf star filling its Roche lobe. The accretion occurs 
along the magnetic field lines down to the magnetic poles 
of the white dwarf. The X-rays and the optical polarized 
fiux, both originating close to the surface of the white 
dwarf are modulated at the orbital period which proves 
the synchronism of the white dwarf rotation with the or- 
bital motion. A comprehensive review of these objects has 
been done by Cropper (1990). 

In the case of BY Cam, the situation concerning the dif- 
ferent periods in the system is not clear. The X-ray source 
H0538-f608 (BY Cam) has been identified as a polar by 
Remillard et al.(1986), on the basis of the detection of 
a circularly polarized optical fiux. This source shows two 
brightness states (Szkody et al. 1990), such a behaviour 
being shared by most polars. However this object is atyp- 
ical by several aspects. In the UV it reveals an abnormal 
emission line spectrum with an enhanced NV line and a 
weak CIV line (Bonnet-Bidaud and Mouchet 1987). This 
could be linked to the chemical composition of a secondary 
whose outer layers have been lost during the evolution of 
the system (Mouchet et al. 1991). This is also reminiscent 
of what is found in some novae, suggesting a possible unno- 
ticed nova- like event (Bonnet-Bidaud and Mouchet 1987). 
However no nova outburst is recorded in archive plates 
(Silber et al. 1992 (hereafter SBIOR)). Noteworthily, X- 
ray spectra of BY Cam obtained with the BBXRT exper- 
iment revealed an oxygen absorption edge near 0.6 keV, 
which intensity either requires an overabundance of oxy- 
gen or partially ionized material (Kallman et al. 1993). 
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Contemporaneous X-ray and optical observations 
(Ishida et al. 1991, SBIOR) have revealed the presence of 
two very close periods around 3.3h, suggesting the asyn- 
chronism of the system (SBIOR). This was the second po- 
lar sharing this peculiarity: the first one being V1500 Cyg, 
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Table 1. Log of the photometric and spectroscopic observations 



Photometry 
date start end phase total exp. 

(UT) (UT) start-end 



V mag. 



Spectroscopy 
start end phase 

(UT) (UT) start-end 



total exp. 



25 Feb. 90 17:40 20:46 0.22-1.15 3h06m 
10 March 91 16:44 21:10 0.67-2.01 3h26m 



14.( 



15-15.6 



12 March 91 16:43 20:13 0.12-1.17 3h30m 14.2-14.9 



17:39 


19:04 


0.22-0.61 


19:51 


20:41 


0.90-1.11 


16:32 


17:42 


0.63-0.93 


17:49 


18:49 


1.02-1.27 


18:57 


20:01 


1.36-1.62 


20:05 


21:06 


1.71-1.96 


16:39 


17:39 


0.11-0.32 


17:44 


18:16 


0.44-0.53 


18:21 


19:33 


0.63-0.92 


19:38 


20:23 


1.01-1.18 



2hl5m 



4hl4m 



3h29m 



a nova which is expected to become synchronous again in 
less than 200 years (Schmidt and Stockman 1991, Katz 
1991). A third one, RX J1940.1-1025, has been added re- 
cently, but contrary to both previous ones, the spin period 
is slightly longer than the orbital period (Fricdrich et al. 
1996). From a compilation of previous optical polarimet- 
ric and photometric measurements of BY Cam, combined 
with new observations, Piirola et al. (1994) obtained a new 
determination of the shortest period (3.3308h) which dif- 
fers significantly from the original determination (Mason 
et al. 1989, hereafter MLS). The presence of two close pe- 
riods in this system led to search for a longer period of the 
order of fourteen days which would be the beat period. A 
large set of photometric data collected over a period of 66 
days, seem to reveal a period of seven days (Silber 1995), 
but a beat period of 14 days has also been suggested by 
Mason et al. (1995a,b). 

Recently, time-resolved UV spectroscopy revealed a 
modulation of the line fluxes and of the radial veloci- 
ties with the longest period identified as the orbital one 
(Zucker et al. 1995), implying an origin far from the ac- 
cretion column, contrary to the common idea that the 
high ionization UV resonance lines are formed close to the 
white dwarf. In the optical emission lines are very complex 
and can show up to four components (MLS). 

In this paper we present an analysis of high temporal 
resolution spectroscopic and photometric optical data ob- 
tained at three different dates. Partial results have been 
presented by Bonnet-Bidaud et al. (1992). The optical 
emission lines, composed of components arising from dif- 
ferent regions, are good tools to constrain the different 
periods present in this system and to trace the geometry 
of the accreting flow. The optical results are compared 
to the UV spectroscopic results. A critical review of the 
different period determinations is also presented in Sec- 



tion 5 and their interpretation is discussed in Section 6 in 
terms of phase-drifts introduced by an asynchronous ro- 
tating magnetosphere. 



2. Observations 

Spectroscopic and photometric data were obtained si- 
multaneously at the AS SAO 6m telescope in Zelenchuk 
(Russia) in February 1990 and March 1991. The obser- 
vations were performed using the SP-124 spectrograph 
and the NEF photometer at the Nasmyth secondary fo- 
cus of the 6 meter Bolshoi Azimuthal Telescope (loan- 
nisiani et al. 1982, Vikuliev et al. 1991). The spectro- 
graph was equipped with a 1200 lines/mm grating. A tele- 
vision scanner with two lines of 1024 channels was used 
to record the sky and source spectra simultaneously in a 
photon-counting mode (see Somova et al. (1982), Drabek 
et al. 1986, Afanasiev et al. 1991, for a detailed descrip- 
tion of the instrumentation). The intrinsic temporal reso- 
lution is 0.033 s. The wavelength range was 3900-4920A, 
4030-5050A and 3950-4970A respectively in Feb. 1990, on 
March 10 and March 12 1991. A wavelength calibration 
has been applied using a He-Ar-Ne lamp. The aperture 
is a circular slit which diameter was chosen between 2 
and 3" , according to the seeing. The corresponding spec- 
tral resolution is '^2.5A. The sky spectrum has been sub- 
tracted but no flux calibration has been applied. The log 
of the observations is reported in Table 1. Photometric 
data were obtained simultaneously with the spectroscopic 
data by splitting the light beam, with about 50% of the 
light being sent into a 12" aperture one-channel fast pho- 
tometer. Continuous Johnson B data were recorded with 
a 0.1s resolution. UBVR were also acquired in one or two 
occasions during each observation and have been used to 
compute the optical magnitude of the source using star F 
from 3C147 field (Neizvestny 1995) as a reference star. The 
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Fig. 2. Enlargements of the optical B light curves of BY Cam 
on March 10 (a) and on March 12 (b) for an illustration of 
the strong variability. The plotted parts are marked A and B 
in Fig.l. The temporal resolution is 1.6 s for both curves. The 
vertical scale is the deviation in % from the mean flux value of 
the corresponding part of the light curve. 



Fig. 1. The airmass and sky-corrected optical B light curves of tion of the different periods of the system, the accuracy 



BY Cam for the three observing dates. The phase is computed 
with the trial ephemeris of Mason et al. (1989). The regions 
marked A and B are enlarged in Figure 2. Note the different 
brightness levels and the strong variability on Feb. 25 1990 and 
March 12 1991. 



given in this ephemeris is purely formal and the period 
(3.322171±0.000017h) is used here only as a trial period 
to facilitate the comparison with other works (see par. 
5.1). 



V magnitude values, reported in Table 1, show the source 
to be in a high state during the three observations, dis- 
tinct from the low state observed by Szkody et al. (1990) in 
January 1989. When two measurements have been done 
during the same night, their different values are consis- 
tent with the amplitude of the light curve (Remillard et 
al. 1986, SBIOR). 

Throughout this paper, heliocentric phases refer to the 
ephemeris given by MLS that is used as a compari- 
son ephemeris. Given the uncertainty in the dctcrmina- 



3. Photometry 

In Figure 1 we show the B light curves obtained simulta- 
neously with the spectroscopic data, at a temporal reso- 
lution of 12.8 s for the three observing dates. They reveal 
two levels of activity. On February 25 1990 and on March 
12 1991 the large amplitude of the modulation is reminis- 
cent of the flaring states observed in the optical as well as 
in the X-rays by Ishida et al. (1991) and SBIOR, while on 
March 10, the light curve shows a lower intensity which 
indicates a level closer to a pulsing state (Ishida et al. 
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1991, SBIOR). For both February 25 and March 12 obser- 
vations, the lowest levels of the light curves which occur at 
the end and at the beginning of the observations respec- 
tively, are consistent with the level detected on March 10. 
The higher fluxes in 25 Feb 1990 and 12 March 1991 com- 
pared to March 10 suggest an additional source of light. 
During a flaring state, SBIOR have reported a onc-pcakcd 
periodic light curve, which exhibits an increase by a factor 
~2 between the minimum and the maximum. The Feb. 
and March 12 curves also exhibit a total variation by a 
factor 2-2.5, but their strong variability excludes a regu- 
lar quasi-sinusoidal shape. However, Silbcr (1995) reports 
flaring state light curves which also depart from a regular 
shape. 

The March 10 light curve varies in amplitude by a lower 
factor of 1.5, quite consistent with the pulsing state light 
curves shown in SBIOR. There is no clear indication of the 
presence of two maxima, while these bumps are usually 
seen in the low-level optical curves with possible unequal 
intensities and separated by a variable phase extension 
(SBIOR, Silber 1995). 

Apart from the large amplitude main modulation, 
these optical light curves show shorter timescale variabil- 
ity. We note a significant (^ 5 niin.) dip seen near phase 
0.1 on March 12. In Figure 2, are reported two blow-ups 
of March 10 and March 12 light curves exhibiting strong 
flaring on a typical timescale of 1-2 minutes on March 10 
and of 3-5 minutes on March 12. The FFT power spec- 
trum analysis of the data does not show significant power 
excess at these periods indicating that these oscillations 
are not coherent. On March 12, quasi-periodic oscillations 
of the order of 30 minutes are visible in the second part 
of the light curve (Fig.l), similar to what reported previ- 
ously by SBIOR. 



4. Spectroscopy 

4-.1. Description of the line profiles 

Spectra averaged over the full observation have been pro- 
duced for the three dates. They are similar to previously 
published spectra, exhibiting the usual strong Balmer and 
helium lines (Remillard et al. 1986, MLS). In order to 
study the orbital and rotational modulation of the line 
profiles, the spectra have been co-added with a 500s reso- 
lution which is the best compromise between a high tem- 
poral resolution and a good signal-to-noise ratio. This re- 
sults in 16, 29 and 22 spectra respectively on 1990 Feb. 25, 
1991 March 10 and March 12. The spectra reveal the pres- 
ence of several variable components in both Balmer and 
He lines. In Fig. 3 the Hell and H^ profiles normalized 
to the continuum are shown for the three epochs in order 
of increasing phase from bottom to top. Large variability 
with phase at a given epoch as well as strong changes be- 



tween the three observations are clearly seen. Remarkably, 
the Hell profile differs from the H^ profile in most spec- 
tra. The most complex profiles appear in Feb. 90 in both 
Balmer and Hell line while in 1991, they are restricted 
to the Balmer components only, the Hell line showing a 
much more regular phase variation. 



4-2. Radial velocity results 

The very complex and variable line profiles observed in BY 
Cam make the radial velocity (RV) analysis very difficult. 
In addition, blends with close weaker lines might occur 
(f.i. Hel 4713A in the red wing of Hell 4686A) and affect 
the measurements. To isolate the different components of 
the complex emission lines, the profiles have been fitted 
with the sum of two or three gaussians of variable widths, 
intensities and positions, using the program SPECTRE 
developed by D. Pclat at Meudon Observatory. In this 
program free parameters are the width, intensity and cen- 
tre of the gaussians as well as the adjacent continuum 
which is fitted with a polynomial of degree one. 

The results of the gaussian fits of the Hell 4686A 
and H^ profiles have been used to derive the velocity 
curves of the different components. The measured radial 
velocities for Hell and for H^ were fitted with a sinusoid 
V{(j)) = J + K sin2'K{(f) — 0o), 4> being the phase com- 
puted from the MLS ephemeris. The best fit parameters 
are reported in Table 2 where the error bars correspond to 
a Icr deviation. The uncertainties given in Table 2 do not 
take into account possible systematic uncertainties in the 
absolute calibration. The errors in the absolute 7-velocity 
are therefore underestimated. 

Multiple components are clearly present and the deconvo- 
lution may not be unique in some cases. The FWHM of 
the different components are given in Table 2, corrected 
for the instrumental response. In February 90, the Hell 
line is wide and exhibits a complex profile: in addition to 
a narrow width (FWHM~200 km s^^) component and a 
broad (FWHM~1140 km s~^) one, a very blue high ve- 
locity shoulder is detectable mostly between phases 0.4 
and 0.6. The RV measurements of the narrow and broad 
components can be satisfactorily fitted with a sinusoidal 
curve, while the measurements of the blue high velocity 
component are more scattered (Fig. 4a,b,c). The narrow 
and broad components appear in phase but with differ- 
ent amplitudes. The high velocity component lags the two 
others by 0.4 in phase. It appears not to be similar to the 
very high velocity component detected by MLS. Indeed 
there is no clear indication of a red counterpart at oppo- 
site phases and the observed modulation is not compati- 
ble with an amplitude as large as r^ 800 km s^^ found by 
MLS. Its low amplitude and its large negative 7- velocity 
seem to favour an origin close to the white dwarf, such as 
a Zeenian component which has been already proposed by 
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Fig. 3. Hell and H^ normalized profiles at the same scale, on February 25 1990, on March 10 1991 and on March 12 1991. 
Phase is increasing from bottom to top as shown on the Y-axis. 



McCarthy et al. (1986) to explain the stationary high ve- 
locity narrow component found in the Hel lines of QQ Vul. 
At the same epoch, the Balmcr lines show more complex 
profiles than the Hell (Fig. 3) and sinusoidal fits cannot 
be quoted unambiguously. 

In March 1991, the Balmer and Hell lines are strongly 
variable. The Hell line essentially exhibits two com- 
ponents, one broad and one of intermediate width 
(FWHM^400 km s^^) with no evidence of the narrower 
(FWHM< 200 km s"^) component seen in Feb. 90. On 
March 10, radial velocities of the intermediate width HeH 
component are satisfactorily described with a sinusoidal 
fit, while on March 12, this component is definitely present 
but its RV curve departs from a sinusoidal curve, being 



quasi-stationary between phase 0.6-1.05. For both dates, 
the radial velocities of the broad component in Hell have 
been fitted, though on March 12 they show a significant 
distortion around phases 0.2-0.4 (see Fig. 4f). The result- 
ing amplitude is lower by 25% on March 12 than on March 
10 and the radial velocities on March 10 lag the ones on 
March 12 by 0.11 in phase (Fig. 4d,f). 
The RV measurements of the broad component in H^ also 
strongly depart from a well-defined sinusoidal curve, par- 
ticularly on March 12 (Fig. 4e,g). The distortions present 
in these RV curves are characterized by a shoulder present 
at phases 0.2-0.5, at the two dates, being markedly pro- 
nounced on March 12. It is reminiscent of what has been 
seen in the flaring state data presented by SBIOR. The 
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Table 2. Averaged deconvolved FWHM and radial velocity sinusoidal parameters of the Hell and H^ line components 



date 


line 


component 


FWHM 
km s^ 


7-ver 
km s^^ 


K amplitude 
km s~^ 


(t>o 


25 Feb. 90 


HcII 


narrow 


200±90 


49±16 


265±24 


0.608±0.013 






broad 


1135± 135 


-50±17 


160±24 


0.591±0.022 






high velocity 


405±70 


-818±15 


158±28 


0.983±0.018 


10 March 91 


H/3 


narrow 


165±35 


18±35 


230±23 


0.138±0.025 




Hell 


intermediate width 


395±100 


-105±10 


172±14 


0.648±0.012 






broad 


1060±260 


-109±10 


357±13 


0.563±0.007 


12 March 91 


H,3 


narrow 


250±45 


-93±51 


214±49 


0.389±0.0027 




Hell 


broad 


1110±220 


-34±10 


250±16 


0.449±0.009 



* These values should be regarded with caution due to the large uncertainty of the absolute wavelength calibration 



RV curves cannot be described by a sinusoidal curve in 
this case. It may indicate either the presence of an addi- 
tional redder component or a partial eclipse of the broad 
component at these phases. Noteworthily, similar devia- 
tions have been reported for VI 500 Cyg by Kaluzny and 
Chlebowski (1988) who have attributed them to the con- 
tribution of a second accreting column. 



nent. Its width and K amplitude are typical of those mea- 
sured in other polars and the origin of this component is 
attributed to regions down to the accretion column (Fer- 
rario et al. 1989, Mouchet 1993b). In 1991 its relative phas- 
ing with respect to the narrow Up component is 0.42 and 
0.06 on March 10 and March 12 respectively, confirming 
the asynchronism of the system. 



At the difference of the Hell line, the Up and H-^ 
profiles show the presence of a narrow (FWHM~150- 
250kms~^) component. The peak is clearly observed 
around phases 0.5 to 0.9 and repeatable from cycle to 
cycle. The velocity curve of this component is shown in 
Figure 5. This component being clearly detectable at spe- 
cific phases only, the corresponding fits are thus less con- 
strained. The amplitude are similar but an obvious lag 
in phase is observed for the two observations separated 
by two days when data are folded with the trial period 
(3.322171h) derived by MLS. 

This phase difference (0.251 ± 0.037) is consistent with an 
orbital period (1.8 ± 0.3%) longer than the trial period. 
The best period derived from our narrow-line observations 
only wiU be (3.3817 ± 0.0085h). This is slightly longer, 
but consistent within the error bars, with the value de- 
rived by SBIOR from the Hq, narrow line velocity. This 
provides an independent confirmation that the period de- 
duced from the narrow-line measurements, usually asso- 
ciated with the orbital period, is significantly longer than 
the period derived from polarimetry. We investigate more 
about the nature of the two periods in Section 5. 

A broad component is obviously present in all data, 
but its characteristics and its corresponding sinusoidal pa- 
rameters can only be safely derived for the Hell compo- 



5. Discussion 

5.1. Emission lines from the horizontal stream 

The true location of the emission line regions in polars 
is still a subject of discussion. Two main components are 
usually considered: a broad line associated with the ac- 
cretion column and linked to the rotation of the white 
dwarf and a narrow line associated with the heated face 
of the secondary, linked to the orbital period. In the case 
of the complex multi-component profiles observed in BY 
Cam, the situation is more intricate and several other 
regions may contribute to build the resulting profile. In 
an effort to reproduce the non-standard observed compo- 
nents, we have investigated a region which has not been 
yet considered in detail, i.e. the horizontal stream of mat- 
ter that links the secondary to the accretion capture re- 
gion, roughly defined by the so-called magnetic capture 
radius. 

Before being captured by the magnetic field lines, 
material forms a stream of matter in the orbital plane 
(Liebert & Stockman 1985) (referred below as the hori- 
zontal stream, as opposed to the accretion column which 
refers to the out-of-plane part of the accreting fiow). This 
stream has been suggested as a possible contributor to 
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Fig. 4. Radial velocities against phase computed using Mason's ephemeris for the different components detected in the line 
profiles. Left: for the three components detected in the Hell profile on February 25 f990: a) broad component; b) narrow 
component; c) blue high velocity component. Middle: for the broad component of the Hell (d) and H/j (e) profiles on March 10 
Right: for the broad component of HeH (f) and H^ (g) on March 12. Best sinusoidal fits are reported (see parameters in Table 
2), except for f and g which have not been fitted since they markedly depart from a sinusoidal curve (see text). 



the emission line profiles (Mukai 1988, Mouchet 1993a). 
The profile characteristics of a component formed in the 
horizontal stream mainly depend on the position and on 
the extension of the emitting region. We have simulated 
the emission from the stream, assuming that its path is 
given by the baUistic trajectory as computed by Lubow 
& Shu (1975). This stream is supposed to extend down 
to a distance Tc at the capture point where the ram pres- 
sure in the horizontal stream equals the magnetic pressure 
(see more details in Bonnet-Bidaud et al. 1996). To evalu- 
ate Tc, we have supposed that the lateral extension of the 
stream is constant, being taken equal to 10^ cm (Lubow 
& Shu 1975, Mukai 1988). We assume a 1 Mq white dwarf 



mass, a typical accretion rate of 10^^^ g s^^, and a dipole 
geometry with a typical polar field of 40 MG consistent 
with the evaluations by Cropper et al.(1989) and Piirola 
et al.(1994). Note however that Mason et al.(1995b) sug- 
gest evidence of a more complex field geometry with ac- 
cretion onto a multipole magnetic field. For the line pro- 
file simulation, free parameters are the mean location of 
the emitting region, its extension and the emissivity of 
each part of the region together with the mass of the 
white dwarf, the inclination angle and the orbital period. 
Since no detailed physical model is available at present to 
give the exact emissivity in the lines, we made the sim- 
ple assumption of an emissivity through the emitting re- 
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Fig. 5. Radial velocities versus phase computed using Mason's 
ephemeris for two observations separated by two days of the 
narrow component of H^ on March 10 (filled circles) and on 
March 12 (open circles). Data are not repeated, March 10 cov- 
ering more than 1.3 cycle. For both dates, the data have been 
offset from the 7- velocity value derived from the fit. The cor- 
responding best sinusoidal fits are shown. Note the significant 
phase lag which can be accounted for by a revised period of 
P=(3.3817±0.0085h). 



gion defined as a gaussian law with a width equal to the 
extension. The emissivity is taken null outside the cho- 
sen portion of the stream. The computed phase-resolved 
profiles have been convolved with a spectral resolution of 
1 A. As illustration, Figure 6 shows different resulting pro- 
files computed at twenty equally spaced orbital phases for 
two selected emitting regions, each defined by their pro- 
jection onto the line of centres and expressed in units of 
the distance X^ between the white dwarf and the inner 
Lagrangian point. Fig. 6a corresponds to a stream exten- 
sion between 0.3 to 0.5 times Xl, and Fig. 6b to a region 
far from the white dwarf, between 0.7 and 0.9 times this 
distance. An inclination angle of 70° has been assumed. 
These profiles clearly show a two-peaked orbitally modu- 
lated width (Fig. 6c) . We stress that the lines produced in 
the stream would therefore in general be strongly variable 
in shape and width, appearing alternatively narrow and 
broad within the same cycle. Such lines would thus be dif- 
ficult to extract by the standard two-component analysis. 
However such high temporal and spectral resolution mea- 
surements are difficult to derive from most data. Therefore 
only the phase-averaged full widths at half maximum ex- 
pected from different zones of the stream have been com- 
pared to the data. In Figure 7 these intrinsic widths are 
reported against the amplitudes K of the RV curves of the 
lines emitted by seven regions of variable positions (a) and 
extensions (b). The inclination angle increases along the 



400 



300 - 



,200 - 




100 - 



Phase 




1.0 


b) 


'JA 






vA 




, 


A 


0.8 


h \ 


rl 




/\l 




lA 


0.6 




J A. 




/a 




/ Aa, . 


0.4 


/ i '■' |l 


^'li 




1 Ml 

;' III 




; 1. 


0.2 


>^\ 


iA 




U. 




h\ 


0.0 


1 \ 





-100(^500 500 1000 
velocity (km s ) 



-1000^500 500 1000 
velocity (km s ) 



Fig. 6. Simulated profiles for emission lines arising in the 'hori- 
zontal' stream, for an inclination angle of 70° and a white dwarf 
mass of IMq. A resolution of lA has been assumed. Left (a): 
extension of the emitting region from 0.3 to 0.5 times the pro- 
jected distance between the white dwarf centre and the inner 
Lagrangian point. Right (b): same as (a) for a smaller exten- 
sion from 0.7 to 0.9 times this distance. Upper panel (c): the 
FWHM versus orbital phase for the two sets of simulated pro- 
files (full line: case a) ; dotted line: case b)). The origin of the 
phase is taken at the inferior conjunction of the secondary. See 
text for more details. 



curves from 10° to 80°. As expected, the wider the emit- 
ting region is, the larger the width is. Maximum values of 
the order of FWHM ~ 400 km s^^, as observed, can only 
be reached for a contribution all along the stream. High 
K values are obtained for the parts of the stream which 
are the closest to the white dwarf. Lines formed in the 
stream at a projected distance from the white dwarf rea- 
sonably larger than 0.1 Xl, the value which corresponds 
to a typical capture radius at ~21 R„,d, have RV ampli- 
tudes restricted to values lower than ^ 800 km s^^ for a 
reasonable (i< 70°) inclination of a non-eclipsing system. 
For higher white dwarf masses (1.4 Mq) the curves are 
shifted up, of the order of 10%, towards greater RV am- 
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plitudes and widths. The plotted curves are computed for 
a system of a 3.33h orbital period, but the values do not 
increase by more than 20% for a typical shorter period of 
2h. 

These stream properties can be compared with the 
characteristics of the components that were identified in 
BY Cam. Narrow components have been found in Balmer 
lines as well as in the Hell line at different epochs (see 
Table 2). The properties of the narrow component mea- 
sured in Feb. 1990 spectra, can be satisfactorily ascribed 
to the stream but its width value is not much constrained. 
A large region starting close to the white dwarf and a low 
inclination angle (i < 30°) can account for the measured 
values. On March 10 1991, the width and the amplitude 
of the H^ narrow line are compatible with an origin in 
the stream, for an inclination angle lower than 30° and a 
large extension of at least O.l-O.SX^ (Fig. 7). The data on 
March 12 are marginally accounted for by an extreme case 
where the complete stream emits. Alternatively, the fact 
that this narrow component is mainly visible at phases 
0.6-1.0 when its radial velocity passes from red-to-blue 
may favour instead an origin in the heated face of the sec- 
ondary (see below). 

More interesting is the intermediate width component 
measured in the Hell line on March 10 1991. Such width 
is not easily produced on the secondary. It can be pro- 
duced in the stream for extended regions and high incli- 
nations, but would have rather high velocity amplitudes. 
The HeH intermediate width component cannot satisfy 
both requirements (K amplitude and FWHM) for an ori- 
gin in the stream, the amplitude being much lower than 
expected from the stream (Fig. 7). In addition, its rela- 
tive phasing with the narrow peak seen in H^ (0.51) is 
far from the expected value if this last component is in- 
deed formed in the hemisphere of the secondary. In fact, 
the simulations show that the blue-to-red phasing of the 
stream component related to the position of the secondary 
is constrained to an interval of phase of 0.75-0.85 depend- 
ing on the exact position, where phase zero refers to the 
inferior conjunction. 

On the basis of these simulations of the line emission in 
the stream, we can also reconsider the possibility that the 
UV lines arise from such a region, as it has been suggested 
by Zucker et al. (1995). The mean K radial velocity am- 
plitude of 370 km s-i for CIV and of 245 km s'^ for Hell 
can be easily reproduced with an origin in the stream. The 
evaluation of the UV line widths is made more difficult in 
the case of the wide NV, CIV and SilV doublets which are 
not resolved with lUE. However, in the case of the Hell 
1640 A line, the blend is made of two very close lines and 
does not affect the determination. Using archive spectra, 
the average measured FWHM of the Hell line is of the 
order of 8 A. The lUE profiles along the dispersion direc- 
tion being well represented by gaussian profiles of about 



5.5±0.3A FWHM at this wavelength (Cassatella et al. 
1985). After deconvolution the intrinsic width of the UV 
Hell line is then of ^ 1060 ± 50 km s~^. The long expo- 
sure time may be responsible for part of the broadening. 
However such a high value is incompatible with the values 
expected from a stream emission (see Fig. 7), unless the 
stream is excessively elongated. We note that due to the 
de-synchronization of the system, the position of the cap- 
ture radius is in fact modulated at the beat period between 
the spin and the orbital periods (see Section 6). We have 
computed this position for given co-latitude angles 9d of 
the magnetic field and for different values of the longitudi- 
nal angle V'd, varying from to 360°, to mimic the variable 
configuration of the magnetic field. We find that the cap- 
ture distance Vc varies very little with 'il)d, the strongest 
modulation (from ^21 to 27 Rmd) being found for a high 
inclined magnetic axis {6d — 85°). We obtain a maximum 
elongated stream for a minimum magnetospheric radius of 
^ 20.6 Rtud, which corresponds to a projected distance on 
the line of the centres of ^0.1 X^. Even an emission from 
a stream extended down to this value cannot account for 
the widths of the lines. An alternative origin of the UV 
lines in the accretion column is proposed in Section 6. 



5.2. The heated hemisphere of the secondary 

As suggested above, the narrow component seen in the H^ 
profile in March 1991 might be thought to arise from the 
X-ray heated hemisphere of the secondary. A lower limit of 
the white dwarf mass can thus be derived using the radial 
velocity amplitude of this narrow component. In the case 
of the extreme hypothesis that the K value represents the 
motion of the centre of the red dwarf, values of the ampli- 
tude K between 200 and 250 km s~^ can be obtained for 
a white dwarf mass> 0.6 Mq and an inclination angle i 
larger than 40°. If, to take into account the decentered po- 
sition of the gravity centre of the illuminated hemisphere, 
we apply a K-correction of 1.3, typical value derived from 
previous studies (Mouchet 1993a), the permitted range of 
WD masses is then 0.8-1.4 Mq, in agreement with the 
lower limit found by Ishida et al. (1991), and i should be 
greater than 50°. 

In contrast, for the UV lines, these large widths exclude 
a dominant origin from the heated secondary. The heated 
hemisphere has been suggested by Zucker et al. (1995) to 
contribute some fraction of the NV line. Indeed a nar- 
row component has been clearly detected in the UV lines 
of AM Her using high resolution spectra (Raymond et 
al. 1995). It has also been invoked to contribute to the 
UV Hell line of V1500 Cyg (Schmidt, Liebert & Stock- 
man 1995), but the formation of this line in V1500 Cyg 
is favoured by the presence of a high temperature white 
dwarf which strongly contributes to the illumination of 
the secondary. 
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Fig. 7. The orbitaUy averaged intrinsic FWHM (in km s~^) versus the amplitude of the radial velocity curves for emission lines 
formed in the 'horizontal stream' for different values of the inclination angle (increasing from 10 to 80° along the curves) for two 
series: at left (a) for a given extension at different positions in the stream; at right (b), for an extension increasing from a region 
close to the white dwarf. The curves are labeled with the corresponding projected extensions of the emitting region, expressed 
in units of X_L,the distance between the white dwarf centre and the inner Lagrangian point. The measurements corresponding 
to the H^ narrow component on March 10 (empty square) and 12 (full square) and to the Hell intermediate width component 
on March 10 (cross) are also reported. 



5.3. Ephemeris: orbital and spin periods 

The determination of an exact ephemeris of BY Cam ap- 
pears to be a tedious task because of the difficulty in 
defining adequate stable time markers in an asynchronous 
binary system. Mason et al. (1989) first determined a 
high accuracy ephemeris based on a recurrent sharp drop 
in the circular polarization data that was assumed to 
track correctly the white dwarf rotation. They noted how- 
ever that their resulting ephemeris was only of appar- 
ent high accuracy, given the highly irregular behaviour 
of the source and the large gaps in their data set. Strictly 
speaking, this first tentative ephemeris with a period of 
(3.322171±1.4 10~^h) can therefore only be used as a com- 
parison ephemeris. A possible alias period near 3.351 h 
was noted by these authors. 

The existence of a different longer period (3.373±0.007h) 
was uncovered by SBIOR from an analysis of the H^ nar- 
row line velocity. A similar value was also derived from 



the line intensity variations. The narrow line emission is 
supposed to track the secondary motion and therefore to 
provide the value of the orbital period. The 1.3% difference 
with Mason's period was then interpreted as an evidence 
of a slight degree of asynchronism of the same type as ob- 
served in the classical nova-polar V1500 Cygni (Stockman 
et al. 1988). 

However as more data accumulate, the difference between 
the two periods becomes less clear. The existence of the 
longer period was confirmed by more data on the H^ nar- 
row component (Sauter 1992, quoted in Zucker et al. 1995) 
and recently by the analysis of the UV emission line veloci- 
ties (Zucker et al. 1995) around a value (3.3558±0.0020 h), 
consistent with Sauter's orbital period. For the shorter pe- 
riod, additional polarimetric data obtained by Piirola et 
al. (1994) were found to be consistent with a short value 
(3.3308196±3.4 IQ-^h), largely outside the error quoted 
in the initial high accurate ephemeris of MLS. A possible 
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Fig. 8. Schematic representation of the different periods found 
in BY Cam with the identification of the origin of the data. 
The open symbols are the different periods determined in this 
paper. Note the distribution of points according to the nature 
of data and the nearly continuous distribution. 



increase of the photomctric/polarimctric period has also 
been suggested by Piirola et al. and Mason et al. (1995a). 
More recently, an even shorter period has been reported 
based on long term photometry by Mason et al. (1995b) 
and Silber (1995) and has been suggested to be a combina- 
tion of the spin and the orbital period. Figure 8 shows the 
different periods presently reported for BY Cam. Ranging 
from 3.29h to 3.38h, they demonstrate the very confusing 
situation for this source. 

In an attempt to clarify this situation, we have per- 
formed a critical re-analysis of existing data, giving partic- 
ular attention to possible alias periods introduced by the 
large gaps between different data sets. These gaps cause 
a loss of the cycle number and prevent the usual (0-C) 
analysis from being performed. We used a minimization 
method in which an estimator is computed from the indi- 
vidual timings as : 
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where P is the trial period, t^, tj the individual tim- 
ings, (Ty the associated errors and N the total number of 
data points. The K factor is introduced to normalize the 
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Fig. 9. a): O-C values of the time arrivals for the radial ve- 
locity measurements of the UV emission lines given by Zucker 
et al. (1995). Dotted line is best linear fit. b): Phase residual 
values after subtraction of the best fit. Note the jump between 
the third and fourth points. The dot-dashed line represents 
the prediction from a phase-drift model for a magnetic dipole 
co-latitude of 12° and a beat period of 14 days (see text). 



estimator to the errors. The best period is defined at the 
estimator minimum. The relevant parameters {tj — ti)jP 
are not independent variables, and thus the estimator does 
not follow a y^ distribution. The statistical significance 
(Prob) was derived from the estimator distribution com- 
puted by a Monte Carlo method. We find that this method 
gives results comparable to the Lomb periodogram (Lomb 
1976) but has the advantage of allowing a weighting of the 
data. 

Using this method, we first tried to combine our opti- 
cal radial velocity narrow line measurements with other 
measurements (MLS, SBIOR). These four collected data 
points span nearly five years in total. The best period is 
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found around 3.38 h, with a 2-day alias at 3.16 h which 
corresponds to the separation of our last two points, but 
a number of aliases close to this best value obviously pre- 
vents a more accurate determination of the period. 
We have also applied this method to the set of data pub- 
lished by Piirola et al. (1994) which gathered timings of 
the circular polarization dips and of a few photometric 
broad minima. In the resulting periodogram, the period 
claimed by Piirola et al. is evident (3.3308±0.0004 h), but 
another longer period at a similar level of significance 
{Prob < 2 10-5) is also present at 3.3749±4 lO"'^ h. We 
note that this last period is consistent with the value de- 
rived by SBIOR. A F-test performed to compare the (0-C) 
residual distributions for the two periods gives no signif- 
icant difference. The polarimetric data collected over ten 
years therefore reveal that both the short and long pe- 
riods coexist in the source. Though at a less significant 
level, a short beat period at a value of 3.2907±0.0004 h is 
also present, compatible with the value derived by Silber 
(1995) and Mason et al. (1995b) from long term photom- 
etry (see Section 6). 

The discovery of a periodicity in the UV line velocity 
at a period of 3.3558 h (Zucker et al. 1995), close to the re- 
fined 3.3554 h long period value of the narrow H^ velocity 
reported by Sauter (1992), was unexpected, since it sug- 
gests an origin far from the white dwarf for the UV lines. 
We have reanalysed the data of Table 1 of Zucker et al. 
(1995). In this case, the accuracy of the data is such that 
one can keep track of the cycles through the full 11 days 
of the observations. A standard (0-C) method can there- 
fore be used instead of a periodogram search and provide 
a more accurate measure of the phase variations inside 
the observations. We have fitted the NV radial velocity 
data of the six separate observations with a sine curve at 
Mason's period (MLS), used as a trial period. The phase 
residuals (0-C) are plotted in Fig. 9. The evident linear 
trend clearly confirms that Mason's period does not cor- 
rectly track the NV velocity. The best period determined 
from the trend is P=(3.3535±0.0029 h), in accordance 
with the value derived by Zucker et al. (1995). However 
after subtraction of this linear fit, significant residual val- 
ues are clearly revealed, defining a saw-tooth shape vari- 
ation through two clearly separated sets of three points. 
The mid-break corresponds to a phase shift of 0.133. If 
the two sets of points had been fitted individually, they 
would have produced two different periods (respectively of 
3.3415±0.0035h and 3.3475±0.0038h for the first and the 
second sets) , each of them shorter than the mean period 
derived using all points. The different "alias" periods are 
shown in Fig. 8, together with all the other period deter- 
minations. The period determinations in BY Cam appear 
roughly distributed in three separate groups depending 
on the nature of the data with long values derived from 
spectroscopy, intermediate values from polarimetry and 
X-rays light curves and short values derived from optical 



photometry. But as illustrated by the UV data, there is an 
apparent continuum of values through the different groups 
(see discussion in Section 6). 



6. Open questions and conclusions 

Though it has been extensively observed, the magnetic 
cataclysmic binary BY Cam is far from been well under- 
stood. The debate about the periods has been reactivated 
by the study of the UV emission lines by Zucker et al. 
(1995). They demonstrate the presence of the so-called or- 
bital period in the radial velocity measurements of these 
lines and are led to assume a formation of the UV lines 
far from the white dwarf, in the orbital plane. We have 
shown above that, because of their broad widths, the bulk 
of these lines cannot be produced either in the horizontal 
stream or in the heated hemisphere. The most natural con- 
tribution to the UV lines is the accreting column out of the 
orbital plane, as it has been previously proposed for sev- 
eral polars, based on the fact that their orbital variations 
are in phase with the broad optical components (Mukai 
et al. 1986, de Martino 1995). We show below that the 
temporal behaviour of the UV lines is indeed consistent 
with this origin. 

The scattering of the periods found in BY Cam (Fig. 8) 
can be explained by a scenario, in which, in the orbital 
frame, the position of the rotational axis is moving slowly 
with the beat period of about fourteen days. The accretion 
column is thus formed along different field lines according 
to the beat phase. Schematically, the accreted material is 
slowly dragged by the magnetic field during the relative 
motion of the white dwarf, up to a certain extent when the 
accretion is no more possible along these lines. Accretion 
has then to occur at the opposite pole, causing the jump 
observed in Fig. 9. 

To compute this effect, we have calculated the position, 
on the white dwarf surface, of the footprints of the field 
lines which intercept the orbital plane at the capture ra- 
dius. This is done for a fixed given co-latitude angle 9d 
of the dipole magnetic field and for different values of the 
longitudinal angle ipd which measure the different phases 
inside the beat cycle. When the capture region is close to 
the white dwarf, the accretion spot is significantly distinct 
from the magnetic pole and for different capture regions, 
its location on the white dwarf surface varies. We find 
that the co-latitude 6^1 of the footprint, with respect to 
the rotational axis, is weakly variable with the beat phase, 
while its longitude ipw, defined in the orbital plane with 
respect to the line joining the centres of the two stars, may 
strongly vary along the beat cycle. This implies a lag of 
the impact spot with respect to the magnetic pole and the 
resulting phase drift would be interpreted as an apparent 
period longer than the true spin period. 
If furthermore, one assumes that the accretion occurs in 
the opposite hemisphere as soon as the threading point 
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is situated at an angular distance from the magnetic axis 
larger than 90°, then at one time, the computed longitude 
abruptly goes down to low values and increases again in 
the cycle. We find that the sudden change of the longitude 
value, due to the switching of the accretion, occurs before 
the longitude reaches the standard value of 180° usually 
expected if one assumes that the accretion occurs at the 
magnetic pole itself. 

In this varying geometry, we have fully computed the 
phasing of the radial velocity curve produced by mate- 
rial falling down the magnetic lines just above the white 
dwarf surface. The amplitude of the shift is determined 
as soon as 8d and the beat period are fixed. The curve 
plotted in Figure 9.b) is computed for values of 9d ~ 12° 
and a beat period of 14 days, and fits the data reasonably 
well. The RV phase slowly varies with the beat period as 
observed, distorting the period determination. The sud- 
den change in the phasing when the pole switches, also 
well reproduces the magnitude of the phase jump {^ 50°) 
observed in the UV 0-C measurements. It appears twice 
during a beat cycle. A beat value can be evaluated from 
the point distribution as 14.5 ±1.5 days. A pole-switching 
behaviour has been also suggested from photometric data 
by Silber (1995, Fig. 3). However large phase uncertainties 
are introduced, in this case, by the fact that the shapes 
of the light curves are variable and strongly depart from 
a sinusoidal curve. 

Strictly speaking, one does not expect to observe the same 
shape of the beat modulation for lines formed close to the 
white dwarf and for lines emitted at other positions along 
the accretion column. The emitting region is also moving, 
depending on the line production mechanism, and may 
introduce additional drifts. The present knowledge of the 
detailed emission line processes at work in these systems 
is not yet sufficient to allow a better evaluation of this 
effect. 

An interesting consequence of this phase-drift model is 
the fact that any period determination is biased depend- 
ing on the length of the observation. Measurements ex- 
tended on more than a beat period will reveal the orbital 
period, while data obtained in a few days will show either 
a shorter period than the orbital one or will not allow any 
period determination if situated close to the jump. Thus 
a large spread of period values may result as it is indeed 
observed in Fig. 8. The considerations above also apply to 
the broad line components usually thought to be formed 
in the accretion column. Interestingly, the radial velocities 
of the Ha broad component measured by Sauter (quoted 
in Zucker et al. 1995) also show variations at the orbital 
period (fi frequency), while they have been found to be 
modulated with the short spin period (lu frequency) by 
SBIOR, based on a set of data spread over six nights only. 
We predict that the 0-C measurements by Sauter would 
mimic the same behaviour as for the NV line. 
In addition our reanalysis of Piirola data has shown that 
a long period is also present in the polarization flux (see 



Section 5.3), together with an indication of the (2lu — Q) 
combination period. These periods are indeed expected for 
a cyclotron emission produced at the basis of the accre- 
tion column (Wynn & King, 1992). By combining the two 
periods of 3.3308h and 3.3749h found in Piirola polariza- 
tion data, a short (2a; — fi) period of 3.2878h is derived 
independently, quite consistent with values determined by 
Silber (1995) and Mason et al. (1995a, b) from photomet- 
ric data. From the same two values, a long beat period 
{lu — 51) of 10.621±3 10^'' days is also derived. This value 
is not strictly consistent with the range of values (13-16 
days) determined from the UV data, and with the similar 
14 day period suggested by Mason et al. (1995a) and Sil- 
ber (1995). 

In conclusion, we have shown that the period determi- 
nations are biased depending on the temporal extension of 
the set of data. The phase-drift model described above ex- 
plains the inability to determine a unique adequate value 
for the periods, when the white dwarf is not exactly syn- 
chronized. This simple picture has to be modified if one 
takes into account more physical complex configurations 
such as a multipole geometry (Mason et al. 1995), a de- 
centered dipole, field line distortions at the threading re- 
gion (Hameury, King & Lasota 1986) or possible inho- 
mogeneous blobs in the infalling material. Moreover it is 
most probable that accretion would occur on both sides 
for intermediate configurations. The orbital period value 
is still inaccurate. It can, in principle, be unambigously 
established from the study of the absorption lines associ- 
ated with the companion atmosphere. A search for the Na 
lines was tentatively done but without success, implying 
that either the companion is very faint or that it is of an 
earlier spectral type than had been expected (Zucker et 
al. 1995). Finally, the discussion of the UV line formation 
suffers from the absence of a 7-velocity value determina- 
tion, which combined with the RV amplitude should allow 
to constrain the emitting region in the accretion column. 
This can be solved with the Hubble Space Telescope and a 
higher spectral resolution than provided by the lUE satel- 
lite. 

Acknowledgments : We are grateful to Didier Pelat for 
providing us his program SPECTRE and to Paul Mason 
for useful discussions. 

References 

Afanasiev V.L., Lipovetsky V.A., Mikhailov V., Nazarov E., 

Shapovalova A.I., 1991, Astrofiz. Issled. (Izv. SAO), 31, 128 
Bonnet-Bidaud J.M., Mouchet M., 1987, A&A 188, 89 
Bonnet-Bidaud J.M., Mouchet M., Somov N.N., Somova T.A., 

1992, in : Stellar Magnetism, Eds Yu.V. Glagolevskij, I.I. 

Romany uk, p 186 
Bonnet-Bidaud J.M., Mouchet M., Somova T.A., Somov N.N., 

1996, A&A 306, 199 
Cassatella A., Barbero J., Benvenuti P., 1985, A&A 144, 335 



14 



The two-period polar BY Cam 



Cropper M., Mason K.O., AUington-Smith J.R., et al., 1989, 
MNRAS 236, 29P 

Cropper M., 1990, Space Sci. Reviews, 54, 195 

Drabek S.V., Kopylov I.M., Somov N.N., Somova T.A., 1986, 
Astrofiz. Issled. (Izv. SAO), 22, 64 

Ferrario L., Wickramasinghc D.T., Tuohy, 1989, ApJ 341, 327 

Friedrich S., Staubert, R., Lamer G. et al., 1996, A&A 306, 
860 

Hameury J.M., King A.R., Lasota J.R, 1986, A&A 162, 71 

Ishida M., Silber A., Bradt H.V. et al. 1991, ApJ 367, 270 

loannisiani B.K., et al., 1982, in : Instrumentation for Astron- 
omy with Large Optical Telescopes, ed. CM. Humphries, 
Reidel, p. 3 

Kallman T.R., Schlegel E.M., Serlemitsos P.J., et al., 1993, 
ApJ 411, 869 

Kaluzny J., Chlebowski T., 1988, ApJ 332, 287 

Katz J.L, 1991, ApJ 374, L59 

Liebert J., Stockman H.S., 1985, in : Cataclysmic Variables and 
Low Mass X-ray Binaries, eds. D.Q. Lamb & J. Patterson, 
Reidel, p 151 

Lomb N.R., 1976, Astrophys. and Space Science, 39, 447 

Lubow S., Shu F., 1975, ApJ 198,383 

de Martino D., 1995, in : Proceedings of the Cape Workshop on 
Magnetic Cataclysmic Variables, eds D.A.H. Buckley and 
B. Warner, A.S.P. Conf. Ser. Vol. 85, p. 238 

McCarthy P., Bowyer S., Clarke J.T., 1986, ApJ 311, 873 

Mason P.A., Liebert J., Schmidt CD., 1989, ApJ 346, 941 
(MLS) 

Mason P.A., Chanmugam C, Andronov LL., et al., 1995a, in : 
Cataclysmic Variables, eds. A.Bianchini, M. Delia Valle & 
M. Orio, Kluwer Academic Publishers, p. 426 

Mason P.A., Andronov I.L., Kolesnikov S.V., Pavlenko E.P., 
Shakovskoy N.M., 1995b, in : Proceedings of the Cape 
Workshop on Magnetic Cataclysmic Variables, eds. D.A.H. 
Buckley and B. Warner, A.S.P. Conf. Ser. Vol. 85, p. 496 

Mouchet M., Bonnet-Bidaud J.M., Hameury J.M., 1991, : 11th 
North American Workshop on CVs and LMXBs, ed. CW. 
Mauche, Cambridge University Press, p 247 

Mouchet M., 1993a, in : Proceedings of the Leiscester Meet- 
ing on White Dwarfs, ed. M.A. Barstow, Kluwer Academic 
Publishers, p 411 

Mouchet M., 1993b, in : Proceedings of the Second Haifa 
Technion Conference: Cataclysmic Variables and Related 
Physics, Annals of the Israel Physical Society, vol. 10, eds. 
O. Regev & C Shaviv, 208 

Mukai K., 1988, MNRAS 232, 175 

Mukai K., Bonnet-Bidaud J.M., Charles P.A., et al. 1986, MN- 
RAS 221, 839 

Neizvestny S.I., 1995, in : Photometric Systems and Standard 
Stars, ed. Moletai Astronomical Observatory, p. 37 

Piirola v., Coyne CV., Takalo S.J.L., Larsson S., Vilhu O., 
1994, A&A 283, 163 

Raymond J.C, Mauche C.W., Bowyer S., Hurwitz M., 1995, 
ApJ 440, 331 

Remillard R.A., Bradt H.V., 
1986, ApJ 302, Lll 

Sauter L., 1992, Ph.D. thesis 

Schmidt CD., Stockman H.S. 



McClintock J.E., Patterson J. 



Silber A.D., 1995, in : Proceedings of the Cape Workshop on 

Magnetic Cataclysmic Variables, eds. D.A.H. Buckley and 

B. Warner, A.S.P. Conf. Ser. Vol. 85, p. 302 
Silber A., Bradt H.V., Ishida M., Ohashi T., RemiUard R.A., 

1992, ApJ 389, 704 (SBIOR) 
Somova T. et al., 1982. in : Instrumentation for Astronomy 

with Large Optical Telescopes, ed. CM. Humphries, Reidel, 

p. 283 
Stockman H.S., Schmidt CD., Lamb D.Q., 1988, ApJ 332, 282 
Vikuliev N., Zinkovskij V., Levitan B., Nazarenko A., 

Neizvestny S., 1991, Astrofiz. Issled. (Izv. SAO) 33, 158 
Szkody P., Downes R.A., Mateo M., 1990, PASP 102, 1310 
Zucker D.B., Raymond J.C, Silber A.D., et al., 1995, ApJ 449, 

310 



1991, ApJ 371, 749 
Schmidt CD., Liebert J., Stockman H.S., 1995, ApJ 441, 414 



This article was processed by the author using Springer- Verlag 
LT^ A&A style file L-AA version 3. 



